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Abstract 



Determination of the top-quark Yukawa couplings in the process e^e^ ttZ 
has been studied for the high luminosity option {J L = 500 fb~ ) of a linear e~^e~ 
collider. Polarization of the electron beam has been considered. The method of 
optimal observables has been adopted to determine the couplings and to disentangle 
different models. It has been found that it could be possible to discriminate between 
the Standard Model scalar Higgs boson, pure pseudoscalar and the Higgs boson of 
mixed CP at about 12a level for ^/s = 1. TeV for the Higgs boson mass slightly 
above twice the top-quark mass both for the leptonic and hadronic top-quark decay 
channels. However, y/s = .5 TeV collider turned out to have too small energy for 
tests of the top-quark Yukawa coupling because of its limited production rate (at 
most 2.5(7 effect). Possible tests of CP violation in the scalar sector have been 
considered. 
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1 Introduction 



It will be very important to directly determine the Yukawa couplings (and there- 
fore CP nature) of any Higgs boson that would be discovered. Although it is very 
plausible that a substantial number of Higgs boson events will first be available 
at the LHC, however future linear e^e~ colliders (NLC) could provide much more 
efficient laboratory to determine couplings of any observed Higgs. The Standard 
Model (SM) predicts an existence of a scalar Higgs boson, however a pseudoscalar 
admixture is absolutely conceivable and as leading to renormalizable interactions, 
it may couple as strongly as the scalar Higgs boson. Since Yukawa couplings 
are proportional to the appropriate fermion masses, obviously the most promising 
coupling to study is the one of the top quark. Here, we are going to discuss a pos- 
sibility to test the top-quark Yukawa coupling in the process e"'"e~ tiZ, our goal 
will be to estimate what would be the future precision for discrimination between 
CP-conserving (scalar and pseudoscalar) and CP-violating Yukawa interactions. 

For our analysis, the relevant part of the interaction Lagrangian takes the fol- 
lowing form: 

C = -!^/ii(a + t^,b)t + c-^^hZ,Z^, (1) 

V 2 cos 0]y 

where g is the SU{2) coupling constant, v is the Higgs boson vacuum expectation 
value (with the normalization adopted here such that v = 2mw/g = 246 GeV) 
a, b and c are parameters which account for possible deviations from the SM 
(a = 1, 6 = and c = 1 reproduce the SM Lagrangian). Since under CP i{a + 
i'~)^h)t ^ i{a — i75&)t and Z^^Z^ ^ Z^Z^ one can observe that terms in the cross 
section proportional to ah or he would indicate CP violation. Since for the process 
considered here {e^e~ — ttZ) there are no ah |^ terms we will focus here on he 
terms. It should be emphasized that this is not just he term but also (6c)^ which 
are signals of CP violation. 

The minimal extension of the SM which provides non-zero h and a, c 7^ 1 is 
the two-Higgs-doublet model |T[] (2HDM). Since the model contains two SU(2) 
Higgs doublets it could be shown [|l| that CP violation may arise in the Higgs 
sector. Such CP violation is explicit if there is no choice of phases such that all 
the Higgs-potential parameters are real. However, even if all potential parameters 
can be chosen to be real, spontaneous CP violation (complex vacuum expectation 
values) is a possibility. Therefore, the 2HDM is a very attractive and simple 
model in which to explore the implications of CP violation in the Higgs sector. 
After SU(2)xU(l) gauge symmetry breaking, if either explicit or spontaneous CP 
violation is present, the three neutral degrees of freedom (one corresponding to the 
Goldstone boson decouples) mix in the mass matrix and as a result of the mixing 
between real and imaginary parts of neutral Higgs fields, the Yukawa interactions 

I* ^ Those terms could appear for polarized ti. 
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of the mass-eigenstates are not invariant under CP. They are given exactly by 
form of the Lagrangian defined in Eq.(|l|). Because of the mixing, the form of the 
ZZh couphng will also be modified by the factor c. Both a, h and c are calculable 
functions of the mixing angles. 

Recently, CP violation in a scalar sector has attracted much more attention as 
it has been realized that it may be generated within the Minimal Supersymmet- 
ric Standard Model (MSSM), where due to the CP- violating soft-supersymmetry- 
breaking Lagrangian complex phases emerge in the effective potential generated at 
the one-loop level of the perturbation expansion Q]. It turns out that the mixing 
between the pseudoscalar A and the heavier Higgs boson H can be as large as 25% 
leading to substantial modification of ZZh\ and Zhih2, where hi^2 are the phys- 
ical mass-eigenstates. Another consequence of the mixing would be CP-violating 
Yukawa interaction of the form given in Eq.(|l]). 

However, in this paper we are not restricting ourselves to any particular model, 
on the contrary, our analysis is supposed to be as model independent as possible 
and we will treat the Lagrangian given by Eq. ([^) as an effective interaction. 

The process e~^e~ —>■ tiZ has been already discussed in the literature by Hagi- 
wara d in the context of the SM. Since the total luminosity for the e+e linear 
collider considered there was a factor of 50 lower than the recently discussed in the 
context of the TESLA project we found it was worth to reconsider the process 
applying some more elaborate statistical methods to test the Yukawa couplings in 
the framework of the effective Lagrangian without referring to any given model. 
Bar-Shalom, Atwood and Soni first noticed in their paper [^] that the process 
e~^e~ tiZ could be used to test CP violation in the Higgs sector, the advantage 
of the process is that CP violation appears even at the tree-level approximation 
of the perturbation expansion. This is why its consequences could be relatively 
large. However Ref. |^ restricts its discussion to the specific case of the 2HDM 
what makes its applicability much more limited. We shall compare our results both 
with Ref. II and Ref. [|. 

The paper is organized as follows. In Section 2 we will briefly review the method 
of the optimal observables. Section 3 will contain a description of our strategy and 
results for the determination of the top-quark Yukawa couplings. Section 4 will be 
devoted to the summary and conclusions concerning tests of the top-quark Yukawa 
couplings at the NLC. In the Appendix, an analytic form for the contribution to 
the matrix element squared linear in the CP-violating couplings will be presented. 

2 The Optimal Analysis Procedure 

We shall consider the process of ttZ production at linear e+e^ colliders: 

e+(p+) + e-(p_) ^ tipt) + iipt) + Z{k), (2) 
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where momenta of the particles involved have been indicated. Since the electron 
beam could be relatively easily polarized, besides unpolarized beam we will discuss 
electrons with positive and negative helicity. The differential cross section for the 
process could be written in the following way: 

= ^ = [fvv + acfac + {acYf(^acf + bcfbc + (fec)V(6c)2] , (3) 

where fyy stands for non-Higgs contribution, whereas other /j are defined through 
coefficients standing in front of them. denotes the final state phase space con- 
figuration, c is the factor which parameterize the ZZh coupling, a and b are the 
scalar- and pseudoscalar-type Yukawa couplings defined in Eq.(p. 

Here, we wish to discriminate between various models on the basis of the di- 
fference between the phase space distribution dependence of fbc, fac, f(bcy and 

/(ac)2 • 

In general, for = J2iCifi{(p) (where are known functions, including 

normalization) the determination of the unknown signal coefficients q with the 
smallest possible statistical error is given ^ by 

Q = ^ Mg'h , where = f ^-^^^l^dcf) , and h = I fk{<P)d<P . (4) 

The covariance matrix for the Cj is 

V,, ^ (AqAc,) = = M^V^eff , (5) 

where cxr = / ^(^0 is the integrated cross section and = Lgsax is the total 
number of events, with L^s being the luminosity times the efficiency. This result is 
the optimal one regardless of the relative magnitudes of the different contributions 
to S(0). It is equivalent to determine c, by maximizing the likelihood (in the 
Gaussian statistics limit) of the fit to the full distribution of all the events. The 
increase in errors due to statistical fluctuations in the presence of background is 
implicit in the possible background contribution to S(0) appearing in Eq. (^), 
which implies larger Mj~^ entries in Eq. (^. From the result of Eq. (^), the in 
the Ci parameter space is then computed as 

XxiY) = E(c: - c^V^\C^ - c,) = E(c: - Q)LeffM,,(c; - c,) , (6) 

where, for our theoretical analyses, the q and Vij are computed within the model 
X, while c- should be calculated for the model Y tested against the model X. The 
value of XxO^) indicates how well could one distinguish X and Y. The power 
of the optimal analysis technique is to take full advantage of all the 
information available in the cross section as a function of the kinematical 
variables. 
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The above procedure is not altered if cuts are imposed on the kinematical phase 
space over which one integrates; one simply restricts all integrals to the accepted 
region. If a subset, 0, of the kinematical variables cannot be determined, then 
the optimal technique can be applied using the variables, 0, that can be observed 
and the functions /j(0) = / fi{<f))d(f). Since cuts restrict an information available 
in a given process, one should expect that the sensitivity of optimal observables 
would be reduced after imposing cuts. 

3 Strategy and Results 

Feynman diagrams describing the process e~^e~ tiZ are shown in Fig. |I[ One can 
think about the process considered here as a way to measure the top-quark Yukawa 
interaction of a Higgs boson produced in a strahlung process off the Z line e^e~ —>■ 
Zh followed by the decay h — ^ ti. However, here we allow also for Higgs boson being 
off-shell, see the diagram (a) in Fig. Within this interpretation the remaining 
(purely SM-type) diagrams (b), (c), (d) and (e) describe a continuum background 
for the Higgs boson decaying to ti pair. Since we are trying to determine the top- 
quark Yukawa coupling the contribution to the matrix element squared relevant 
for us would be the interference of the Higgs-boson-exchange diagram (a) with Z- 
boson radiation off electron line (b), (c) and top-quark line (d), (e) and the square 
of the Higgs diagram (a). In some sense the process e~^e~ tiZ considered here 
is complementary to e^e^ tih discussed recently in Ref.[0, as both processes 
originate from the Higgs-boson strahlung off Z line followed either hj h —>■ ti or 
Z ti, respectively. 

Hereafter we will be using the top-quark mass rrit = 175 GeV, Z-boson mass 
niz = 91.17 GeV and the Weinberg angle sin^ 9 = 0.23. Electrons will be considered 
as massless. Before we start the optimal analysis it is useful to recollect the SM 
predictions for the total Higgs-boson width and also the total cross section for 
the process of interest {e~^e~ — >■ tiZ) for various ^/s and Higgs-boson masses M^, for 
polarized j = ±1 and unpolarized j = electron beams. Results are presented in 
Fig. 1^0. Because the Higgs boson could be on its mass shell, the Breit-Wigner form 
of the propagator must be adopted, therefore it is relevant to remind ourself the 
strong increase of Th as a function of Mh, see plot (a) in Fig. |^. Concerning the total 
cross section as a function of ^/s we observe a threshold increase at ^/s = mz + ^rrit 
followed by a typical s-channel suppression for higher energies. For the polarization 
dependence, the relative suppression of the polarization j = I could have been 
anticipated, as it is caused by a constructive and destructive interference between 
7 and Z exchange diagrams for j = —1 and j = +1, respectively. As it is seen from 

our code we have obtained exactly the results presented for the SM by Hagiwara et al. 
in Ref . . Other test we have performed was to reproduce the narrow- width-approximation for 
the on-shell Higgs boson. We have also checked our code against the analytical form shown in 
the Appendix for the contribution linear in be. 
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plots (c) and (d) in Fig. |^ the relative contribution from the Higgs boson to the total 
cross section could reach up to 30% in the region where 2mt < < \fs — mz- It 
could be read from Fig. |^ that for a given ^/s the total cross section for low Higgs- 
boson mass, = 200, is greater than the one for high mass. That means that 
the Higgs-boson-exchange diagram (a) adds constructively with the other diagrams 
for low masses and interferes destructively in the region of high masses. In fact, 
it is seen from Fig. ^ where separate contributions to crtot(e"'"e~ — > tiZ) have been 
plotted. This is related to the correct high-energy behavior of the cross section 
in the large limit. Since the virtual Z boson couples to massless electrons, 
effectively it is mostly the transverse component which propagates and therefore 
the Higgs-boson exchange is not necessary to provide "good" high-energy behavior 
in the infinite Higgs-boson mass limit. 0. Quantitatively however, the effect is 
totally negligible for the total cross section. 




Figure 1: The Feynman diagrams for the process e~^e ttZ. 



Following detector design of the Japanese Linear Collider |^ we impose here 
the polar angle cut | cos 6*1 < 0.98. We will consider separately two cases. For 
the first one, which we call the leptonic case, one top quark decays hadronically 
and the other leptonically and the second one, which we call the hadronic case 
in which both t and i decay hadronically. We choose the following ti tagging 
efficiencies: .30 for leptonic and .44 for hadronic modes. For hadronic decays it 

"^For a detailed discussion of unitarity constraints see Ref. O]. 
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is very difficult to disentangle experimentally t and i, we would be only able to 
tell from the invariant mass reconstruction that a given jet corresponds to t or 
t. Therefore, adopting the optimal observables for hadronic ti decays, we will 
always be averaging over t and i that would correspond to t and i tagging just 
by the invariant mass measurement without determining the charge of decaying 
top quarks. Obviously, for leptonic decays it is not necessary since the charge of 
hadronically decaying top quark would be identified by the lepton charge from 
the semileptonic decay of the other top quark. For simplicity we will assume 100% 
efficiency for Z tagging. For definiteness, in what follows, we will consider one year 
of operation for high luminosity option of the NLC that has been examined in the 

V 



context of the TESLA collider design, for which one expects L = 500 fb ^ 
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Figure 2: The SM predictions for the total Higgs-boson width (a), the total cross 
section for e^e~ — ttZ for indicated polarization of the initial electron beam 
[j = — 1,0, +1) as a function of ^/s for fixed Higgs-boson mass (b), and the 
total cross section as a function of Mh for y/s = .5 TeV (c) and for -y/s = 1. TeV 
(d). 
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Figure 3: The SM predictions for the total cross section for e'^e~ — > tiZ for 
indicated polarization of the initial electron beam (j = —1, 0, +1) as a function of 
Higgs-boson mass for ^/s = .5 TeV and for ^/s = 1. TeV. Contributions from 
the Higgs-boson diagram squared (hh), all non-Higgs diagrams (VV), and from 
the interference between Higgs and non-Higgs diagrams (hV) have been plotted 
separately. The interference becomes negative for Mh > 390 GeV and 700 GeV 
at y/s = .5 TeV and = l.TeV, respectively. Therefore we have plotted the 
absolute value of the interference term. The solid line stands for the sum of all the 
contributions. 

We first focus on the detection of CP violation in the process e^e^ — > ttZ. As 
has already been mentioned CP, violation would result in the appearance of non- 
zero be and {bcY terms in the differential distribution, Eq.(^). An interference of 
the Higgs-exchange Feynman diagram with the Z radiation off the final top-quark 
lines generates terms be, whereas (6c) ^ are emerging from the square of the Higgs- 
exchange diagram. Let us assume that the non-standard CP-odd admixture to the 
Yukawa coupling is small, so that the terms proportional to (6c)^ could be neglected 
and the leading contribution is the one linear in be. In the Appendix we present 
the analytical formula for the contributions to the matrix element squared that 
are proportional to be. It turns out that in the small be limit the differential cross 
section could be written in the form S(0) = Ilieifi{(f)) and the optimal analysis 
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outlined in the previous section could be directly applicable 0. 
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Figure 4: The statistical error 6{bc) for the determination of be for the case when 
{bcY could be neglected, calculated for the indicated polarizations of the initial 
electron beam and ^/s = .5, 1. TeV. 



As it has been shown in Ref. one can construct an appropriate weighting 
function uji{(f)) such that / uJi{(f))Tj{(j)) = Ci. The covariance matrix for such ob- 
servables is given by the Eq.(^. In our simple case when only linear terms in the 
unknown coupling be are kept, the weighting function for determination of be and 

"''Notice that the Higgs-boson width Th enters the matrix element squared through the de- 
nominator of the Higgs boson propagator. Since Th depends on a^, 6^ and the differential 
cross section is not a polynomial in couplings to be determined, therefore in general Eqs.(^,^) are 
not applicable. However, for small b one can expand the propagator around its SIVI form. Since 
we are going to keep only linear terms in 6, effectively we are allowed to use the SIVI form of 
the propagator while calculating the functions fi. Therefore, in the small b limit, the differential 
cross section is of the desired form indeed. 
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its statistical error become: 



/ 




(7) 



1 1/2 



5{bc) 



1 / ^SM 

N r fl 



(8) 



where T^sm denotes SM contribution to the differential cross section. Numerical 
results for 5{hc) are shown in Fig. ^ As seen from the figure, the smallest errors 
appear for such that the Higgs boson could be on its mass shell and decay to 
tt, this is the region corresponding to the maximal cross section, see Fig. In 
general, hadronic top-quark decay modes provide much larger errors, however in 
the vicinity of = 370 GeV both hadronic and leptonic final states lead to errors 
of the order of 1. for ^/s = .5TeV. For ^/s = 1. TeV minimal errors are obtained 
for Mh = 450 GeV for leptonic decays and one can see that leptonic modes offer 
S{bc) ~ .25 whereas for the hadronic ones we get S{bc) > 1.. The enhancement of 
the error for the hadronic modes illustrates how important the t — t identification 
is. As expected, loosing information on the top-quark charge leads to larger errors. 

Let us now consider the general case with non-negligible terms 0{b'^). Then the 
dependence of the differential distribution Eq.(^) on the unknown couplings a, b 
and c is much more involved since the Higgs-boson width, which enters the denom- 
inator of the Higgs-boson propagator also depends on those couplings. In other 
words, functions fi appearing in the distribution shown in Eq.(|^) are unknown as 
being a, b and c dependent. However, we shall assume here that at the time when 
the NLC will operate the total Higgs-boson width will be known from mea- 
surements at the LHC Therefore we are going to consider only those extensions 
to the SM that predict the same F/j, otherwise the model would be experimentally 
excluded at the very beginning. Moreover, the measurement of the total cross 
section for the process e~^e~ ttZ would also eliminate those models for which 
the cross section would differ from the measured one. Therefore, in addition, we 
will restrict ourself to those models which provide the total cross section identical 
with the SM one [5 The models we are discussing are specified by the couplings 
a, b and c appearing in the Lagrangian given by Eq.(|l|), the following versions will 
be considered: 

• SM: a = 1, 6 = 0, c = 1 

• CP-mix. 1: a = b > 

• CP-mix.2: a = b<0 

• CP-odd: a = 
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Figure 5: Solutions for a and 6 as a function the the Higgs-boson mass. For a given 
Mh, plotted a and b reproduce SM values for Th and crtot(e"'"e~ — * ttZ) for various 
polarizations of the initial electron beam. 



For given Mh and values, the width Th and the cross section (Ttot(e^e^ 
ttZ) could be read directly from Fig. 0. Specific values for a, b and c for the last 
three models are obtained by solving the quadratic equations resulting from our 
requirements of Th = Ff^'^ and atotie^e~ tiZ) = af„f (e+e" tiZ) 0. As the 
equations to be solved are quadratic we have several solutions, we have chosen to 
discuss only those with the value of c closest to 1. In the CP-mixed case there are 
both solutions with a = 6 > and a = 6 < we chose to discuss both cases. Since 

''^The choice of the SIVI values for the width and the cross section is, of course, arbitrary. It 
serves just a possible example of the future experimental data. 

"^It should be mentioned that relaxing this constraint may lead to different conclusions, see 
Ref. j^, where models corresponding to very distinct and (T(ot(e~''e^ tiZ) have been com- 
pared. We have checked that for models compared in Ref. [|| either widths or total cross sections 
differ by about a factor of 2. 
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for all considered models it turns out that c is very close to 1 (for Mh < \fs — mz we 
have found 0.98 < c < 1.05, for higher masses c = 1±.07) we present in Fig. ^only 
solutions for a and 6 as a function of Mh. As it is seen, substantial deviations from 
the SM are possible keeping Yh = T^^^ and atotie+e' ttZ) = a^^f (e+e" ttZ). 
Since the total cross section (see Fig. ^ depends on the initial beam polarization 
the solutions for a and b also vary with the polarization. Kinks in the curves 
correspond to the region where Higgs-boson contribution to the total cross section 
for (Tto<(e^e~ ttZ) is maximal, i.e. in the vicinity of the threshold for a real 
Higgs boson decaying into ti. 
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Figure 6: Xx(^) = yXxO^) where XxO^) defined by Eq.(|D calculated for 
semileptonic top-quark decay final states (at least one W boson decays lep- 
tonically). The models considered are specified in the text. Results are pre- 
sented for polarized initial electron beams (j = —1,0,-1-1) and the total energy 
v/i= .5, l.TeV. 
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Since we stay within a class of models predicting the same |^ Th, 
therefore we are entitled to adopt the method of the optimal observables 
since the functions /j from Eq.(^) are known indeed. In Figs. |^ and ^ we 

show results for XxiY) = \/x'xO^) calculated according to Eq.(|^) for the leptonic 
and hadronic events, respectively. First of all we can see that for ^/s = .5 TeV even 
if the real Higgs boson is decaying into tt, x never exceeds 2.6 both for leptonic 
and hadronic top-quark decays, making any conclusions very uncertain. The initial 
electron polarization may be helpful to some extend, although its effect is never 
very large. For i/s = 1. TeV we observe for the leptonic events that xsMiCP — odd) 
is reaching 12.0 just above mu = ^nit and it is above 3.8 up to Mh = 600 GeV. For 
the hadronic events the maximal x is more or less the same, however the region 
of large x is narrower. This should be understood as an consequence of averaging 
over t and t as it suppresses contributions from terms linear be leaving just the 
contribution from the square of the Higgs-exchange. 

It is worth mentioning that even though Fig. |^ was generated for the SM it 
is, in part, still applicable in the context of the extensions defined above. This is 
because the constraint Th = implies for c ~ 1 that Th{h tt) should also 
be of the SM size. Therefore, approximately, the contribution from the Higgs- 
boson diagram squared should be close to the SM one for Mh > 2mt Q Since 
c"iot(e'''e~ — > tiZ) = (yf^{e^e~ — >• tiZ), it is obvious that the interference between 
the Higgs and non-Higgs diagrams must also be the same. 

One could notice that, in spite of the fact that atotU = +1) < o'totij = 0) < 
c"tot(j = ~1) < there are such models that for certain Higgs-boson masses we have 
x{j = +1) > x{j = 0) > x{j = Oils should however remember that for a 

given Mh the model parameters a, b and c adopted for j = ±, are not the same as 
they are adjusted to reproduce the SM prediction for Th and cTtot(e+e~ tiZ), see 
Fig^. The functions fi utilized according to Eq.(^), for the calculation of Mij are 
also polarization dependent. In consequence it may happen that the polarization 
leading to a lower total number of events provides higher x for certain models. 

One could also try to disentangle the CP-odd model and the models of mixed- 
CP- Yukawa couplings. However it turns out that the Xcp-oddiCP — mix.l) is 
always below 0.4 and 2.3 for -y/i = .5 and 1. TeV, respectively. In other words 
the models are too similar to be distinguished even by the method of optimal 
observables. As for Xcp-odd{CP — mix.2) the maximal values obtained at y/s = 
.5 and 1. TeV are 1.0 and 6.2, respectively. Leptonic decay modes led always to 
higher x while distinguishing the CP-odd model and models of mixed-CP-Yukawa 
couplings. 

It should be mentioned here that if rriH turned out to be just above 2mt then 
non-zero top-quark-width effects would smear out the "spikes" observed in Figs. ^, 

"^The one measured in the future at LHC. Here it is assumed to be F^*^. 
"^As it is saturated for Mh > 2mt by the Higgs-boson resonance contribution : tTtot(e+e^ 
ttZ) ~ atot{hZ)Th{h ^ tt)/Th 
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Figure 7: XxO^) as in Fig. |^ but for purely hadronic top-quark decays. 

0, by an amount or the order of the width ~ 1.5 GeV. In that case one can 
expect the QCD corrections to be relevant, of the order of 25% Q. 

So far we have been considering only production of unpolarized tt pairs. How- 
ever, it is well known that helicity of the top quark could be statistically determined 
through its semileptonic decay modes, e.g. lepton energy distribution is sensitive 
to the initial top-quark helicity ||^. It is easy to notice that for polarized ti pairs 
differential distributions would contain, in addition to be, also other possible signal 
of CP violation, namely terms proportional to ab. To perform a realistic analysis 
one should consider t, i and Z decays and then having 8-body final phase space 
determine ab and be applying the method of optimal observables, this is however 
beyond the scope of this paper. 

''^For the analogous process e+e" tih the QCD corrections are known and for mu = 
100 GeV ~ mz they have been found to be at the level of 35%. 
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4 Summary and Conclusions 



The method of the optimal observables @] has been apphed to study Yukawa 
couphngs of the top quark in the process e~^e~ — tiZ. The high luminosity 
(L = 500 ih~^y~^) option of linear e~^e~ colliders examined in the context of 
the TESLA collider design has been adopted. Using the optimal weighting 
function uJbc given by Eq.(^ we have found that the minimal statistical error for 
establishing CP violation in the scalar sector through determination of be, 6 (be) 
could be of the order of .25 for Mh — 450 GeV for the ^/s = 1. TeV collider. 
The optimal observables have been used to disentangle the SM (CP-even) and 
CP-odd or CP-mixed Yukawa couplings. In order to make analysis more realistic 
only models predicting the total Higgs-boson width and the total cross section for 
e~^e~ — > tiZ equal to those of the SM have been considered. It turned out that 
the lower energy option ^/s = .5 TeV of the collider provides too small production 
rate to be applicable for the determination of the top-quark Yukawa coupling. For 
^/s = l.TeV one could distinguish the models considered here even at the level 
of X = 12 (for Mh — 350 GeV). However, the results presented in this paper 
correspond to nearly ideal experimental setup. The full analysis should include 
background and detector simulation which would not only lower the efficiencies 
but also, due to the finite detector resolution, would smear both the signal and 
background distributions. A detailed study including all those effects would be 
desired. Since the statistical significances found here were relatively large one can 
believe that in a real experiment it would possible to determine the top-quark 
Yukawa coupling in the region where 2mt < mu < 600 GeV for y/s = 1 TeV. 
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Appendix 

Below we present contributions to the matrix element squared which are linear 
in be. As it is seen we have separated terms proportional to the real (ReH/j) and 
imaginary (ImH/j) part of the Higgs-boson propagator Ilh{pD = \p\—m\+imhT 
for ph = pt+ pf. 

(1 7D 
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where 
Iv 



+9ht9vum - fh)m% + ~g{p - {h - sf + m|s)] ] + 

e\h'^ + {2h - s){ml - s)] - fg(h + 2m\] 
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It was useful to adopt above the following CP-even and CP-odd quantities: 

e = {p_- p+) ■ {pt - pi) 

f = {P-- P+) ■ {Pt + Pt) 
g = (p_ + p+) ■ {pt - Pt) 
h = {p-+ p+) ■ {pt + Pt) 



(10) 



The vector-boson propagators are defined as ny(s) = [s 
propagators the following notation has been used: 
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where Uxiq"^) = [q 
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Table 1: The vector and axial- vector couplings of the electron and of the top quark. 



The contribution proportional to the real part of the Higgs propagator have 
been already published in the literature ^ . We have noticed two missprints present 
in the published result 0: factor 1/2 is missing in Eq.(12) and sign in front of the 
first term in Eq.(13) should be reversed. 

It has been argued in Ref. that the contribution from the imaginary part 
of the Higgs propagator should be omitted since it is of higher order as being 
proportional to the Higgs width T^. However, one should have in mind that the 
integration over the phase space around the Higgs-boson mass provides an extra 
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factor which cancels Th from the numerator of the imaginary part of the prop- 
agator leading to the contribution of the same order as the one proportional to 
the real part. This is exactly the same mechanism which is responsible for the 
applicability of the narrow width approximation. However, to generate the com- 
pensating factor it is necessary to have the pole of the Higgs-boson propagator 
with its vicinity (of the size of the Th) within the integration region. In the process 
under consideration, the square of the Higgs boson momentum, pi = {pt + Pt)^ is 
limited by 4m^ and s — m|. It is easy to find that the position in the middle of 
the allowed region corresponds to rrih = 430 and 750 GeV for = .5 and 1 TeV, 
respectively. Since both for 430 GeV and 750 GeV the allowed region contains the 
Higgs-boson pole with its vicinity of the size of the width one can expect contri- 
butions of the same order from both real and imaginary parts of the Higgs-boson 
propagator. We have found that the effect of the imaginary part could increase 
Xcp-odd{C P — mix.l) even by 200%, however for that case x was still below 1. The 
only phenomenologically relevant case for which the imaginary part was important 
was xsm{C P —mix.2) for — 500 GeV where we have found 10% increase caused 
by the imaginary part. It turns out that for all other considered cases the contri- 
bution from the imaginary part of the Higgs-boson propagator could be neglected 
for all practical purposes, although a consistent treatment demands both the real 
and imaginary contributions. 
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